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Abstract 

Accurate calculation of the motion of a compact object in a background spacetime induced 
by a supermassive black hole is required for the future detection of such binary systems by the 
gravitational- wave detector LISA. Reaching the desired accuracy requires calculation of the second- 
order gravitational perturbations produced by the compact object. At the point particle limit the 
second-order gravitational perturbation equations turn out to have highly singular source terms, for 
which the standard retarded solutions diverge. Here we study a simplified scalar toy-model in which 
a point particle induces a nonlinear scalar field in a given curved spacetime. The corresponding 
second-order scalar perturbation equation in this model is found to have a similar singular source 
term, and therefore its standard retarded solutions diverge. We develop a regularization method 
for constructing well-defined causal solutions for this equation. Notably these solutions differ from 
the standard retarded solutions, which are ill-defined in this case. 
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I. INTRODUCTION 



One of the interesting sources of gravitational waves which could be detected by the 
planned laser interferometer space antenna (LISA) is composed of a compact stellar 
object with a mass yU (e.g., a neutron star or a small black hole) which inspirals towards 
a supermassive black hole with mass M, such that M ^ fi. Successful detection of these 
sources using matched-filtering data-analysis techniques requires preparation of gravitational 
waveform templates for the expected gravitational waves. The calculation of a complete 
waveform template corresponding to an inspiral orbit which lasts for one year, requires 
calculation of about 10^ wave cycles P]. One of the main difficulties in the construction of 
templates for such long wave-trains is to keep track of the gravitational-wave accumulating 
phase. Accurate calculation of this phase requires accurate calculation of the compact 
object inspiral trajectory, which follows from the solution to Einstein field's equations for 
this problem. 

In practice one can take advantage of the smallness of the ratio fi/M and use perturbation 
analysis to simplify the calculation of the inspiral trajectory. For concreteness suppose that 
the small compact object is a small Schwarzschild black hole, which moves in a vacuum 
background spacetime induced by a supermassive black hole. Here, far from the small black 
hole (i.e., at distances much larger then fi) the geometry is dominated by the supermassive 
black hole. In this region it is useful to express the full spacetime metric g^^ as a sum 
of a background metric Qf^^, induced by the supermassive black hole, and a sequence of 
gravitational perturbations that are induced by the small black hole, reading 

g,Ax) = 9,uix) + fih^^^ix) + f^XK^) + Oif^') . (1) 

Here the perturbations {h^fll} are defined on the background geometry, and are independent 
of fi. By substituting this asymptotic expansion into Einstein's field equations one obtains 
linear field equations for the first-order and second-order gravitational perturbations /i^^ 

(2) 

and hfj,,), respectively. 

To further simplify the calculation it is useful to consider the limit /x — of the series 
(H)). Notice that by definition h^/^u and /i^S do not depend on jj and therefore the form 
of the perturbations equations which they satisfy is not affected by this limit. However, 
the domain of validity of these equation is in fact expanded as /i — 0. For a finite value 
of fj, the domain of validity of expansion (jT]) hes outside a worldtube of radius rE(/i) such 
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that te ^ fi- At the hmit /i — ^ the radius te approaches zero and the series ((T)) is 
vahd throughout the entire background spacetime excluding a worldhne 7. Here D'Eath has 
shown that in this hmit the first-order metric perturbations h^f}} (in the Lorenz gauge) 
are identical to the corresponding perturbations that are produced by a unit-mass point 
particle tracing the same worldline 7. Therefore, for practical purposes (e.g., to calculate 
the emitted gravitational waves) we can identify 7 as the worldline of the small black hole in 
the limit fi ^ 0. Note that for a finite value of fi this identification is valid only at distances 
from the small black hole which are much larger then /i. 

Complementing the metric perturbations series, we also treat the worldline 7 perturba- 
tively. Here at the leading order of approximation the worldline 7 is dominated by the 
background geometry. In fact it was found to be a geodesic in the background spacetime 
(see e.g., |3|). 

At the next order, the interaction between the small black hole's mass and the grav- 
itational field induced by the small black hole itself, gives rise to a gravitational self-force 
acting on the small black hole. The leading order effect of the self-force originates from 
interaction with the first-order gravitational perturbations This leading self- force (of 

order yU^) produces an acceleration of order fi for the black hole's worldline 7. For a vacuum 
background geometry, formal and general expression for the leading order self-force on a 
small compact object was derived by Mino, Sasaki, and Tanaka Q], and independently by 
Quinn and Wald ^ using a different method^ Later, practical methods^ calculate this 



■jater, practical metnoas to 
1 3,0 U, l3; see also Q 



self- force were developed by several authors p, |7|, |8|, |9|, |10|, lUl, 1121 ; see also [13[ for a different 
approach to this problem. 

The next order correction to 7 in this approximation scheme originates from the interac- 
tion of the black hole's mass with the second-order gravitational perturbations fi'^h^^^. Here 
it was shown by Poisson Q| by simplifying an argument by Burko that the shift in 
the inspiral orbit due to this interaction may give rise to a shift in the gravitational-wave 
phase of order This 0(1) shift in the phase is important for the construction of long and 
accurate waveform templates for LISA, and provides a practical motivation for the study of 
the second-order gravitational perturbations. 

Schematically (and without indices) the second-order perturbation equations, which are 
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obtained from Einstein's field equations in the limit /i — 0, read 



(2) 



Here D denotes a linear differential operator, and S denotes the corresponding source term. 
Recall that in the limit of interest we may consider h^-^^ as being induced by a unit mass 
point-particle. Therefore, in the Lorenz gauge this field satisfies a linear wave equation 
with a point-particle source term. The solution of this wave equation has a well known e^^ 
divergent behavior in the vicinity of the worldline 7, where e denotes the spatial distance 
from 7. (Throughout this paper we shall define the spatial distance e of a point x from 
a worldline to be the length of a spacelike geodesic which connects the worldline with x. 
This geodesic is defined such that its tangent vector is normal to the worldline, and it lies 
within a normal neighborhood of this worldline). Examining the form of the source term S 
reveals that it contains two types of terms that can be schematically expressed as V/i*^^^ V/i*^^^ 
and /i*^^) VV/i^^-*, where V schematically denotes the covariant derivative with respect to the 
background metric. Here there is a problem, since the 0(e^^) divergent behavior of h^^^ 
implies that source term S diverges like e~'^. A naive attempt to construct a standard 
retarded solution to Eq. say by imposing Lorenz gauge conditions on h^'^^ and then 

formally integrating S with the corresponding retarded Green's function, ends up with an 
ill-defined integral for the retarded solution which diverges at every point in spacetime. 

In this article we lay the foundations for a regularization method which can be generalized 
to enable one to construct a well-defined solution to Eq. (0)- The construction of a well- 
defined solution is a first step towards the calculation of the previously mentioned 0(1) 
corrections to the gravitational- wave phase-shift, which are produced by the interaction 
between the second-order gravitational perturbations and the small black hole. 

In this paper we study a simpler nonlinear scalar toy-model, which exhibit similar diffi- 
culties. We consider a point-particle with a scalar charge q which induces a scalar field 
that satisfies a certain nonlinear field equation in a given curved background spacetime. Re- 
call that in the corresponding gravitational problem we found it useful to consider the limit 
/i — > 0. Similarly in our scalar model we shall seek a solution to the nonlinear scalar field 
equation in the limit q ^ 0. For small values of q we shall derive linear partial-differential- 
equations for the first-order and second-order scalar perturbations. In the limit the source 
term in the second-order scalar perturbation equation turns out to be highly singular, and 
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diverges like e~^. This singularity is similar to the singularity encountered in the source 
term of the corresponding second-order gravitational perturbations equation Similar to 
the gravitational case, the standard retarded solution of the second-order scalar perturbation 
equation diverges at every point in spacetime. To resolve this difficulty, we develop a new 
regularization method. Essentially this method is based on constructing a new well-defined 
causal solution to the problematic second-order scalar equation. Notably this solution is 
different from the standard retarded solution which is ill-defined for this problem. The 
generalization of the regularization method presented here to the more physical, and more 
difficult second-order gravitational perturbations will be published elsewhere. 

This article is organized as follows: First, in Sec. Owe present our nonlinear scalar toy- 
model, and construct the relevant perturbations equations; then in Sec. Illll we construct 
the first-order scalar perturbation in the limit q —>■ 0; and finally in Sec. IIVI we consider 
the second-order scalar perturbation equation, for which we construct well-defined causal 
solution in the limit g ^ 0. 

II. NONLINEAR SCALAR TOY-MODEL 

Consider a particle with a scalar charge q which traces a geodesic ^ z{t) in a given 
background spacetime. Here r denotes the particle's proper time, with respect to the back- 
ground geometry, the background metric is denoted with gf^,^, and we adopt the signature 
(— , -|-, -|-, -|-) throughout. In our model, the scalar field satisfies the following nonlinear 
field equation 

□0 = -Anqn + V^^V^ . (3) 

Here □ denotes the linear wave operator given by □</> = g'^'^V pcj), and n{x) denotes the 
particle's number density, reading 

/"°° 1 

n{x)= ^=5^[x - z{T)]dT , (4) 

J-oo V-9 

^ Note that the problematic 0{e^^) source term in the corresponding gravitational equation (0) originates 
from terms involving h^^\ where h^^^ is induced by a unit mass which traces a geodesic worldline. This 
motivates choosing a geodesic worldline for the scalar model as well. To keep our scalar particle on a 
geodesic trajectory we assume that external forces are acting on the particle such that the sum of all 
the forces which act on the particle vanishes - this sum of forces includes the scalar self-force and all the 
external forces. 



5 



where g denotes the determinant of the background metric, and S'^[x — z{t)] formally rep- 
resents a distribution which is similar to the four- dimensional Dirac delta-function in the 
sense that it vanishes away from the worldline. While the Dirac delta-function is an ex- 
tremely useful mathematical tool to represent a structureless point particle in the context 
of a linear field equation, we find that it is somewhat less suitable for the nonlinear field 
equation discussed in this paper. An alternative mathematical formulation which precisely 
encodes the physical properties of the source is provided below. 
We seek a perturbative solution to Eq. (jS)) in the form of 

In this expansion the perturbations (f)^^^ , 0*^^^ are independent of q. Since is a dimensionless 
field, the charge q has to form dimensionless combinations with other dimensional quantities 
in each term in Eq. (0). The relevant dimensionless combinations are q/e and {q/TZi}, where 
{TZi} denotes the set of length scales which characterize the background spacetime. Note 
that expansion is valid only if these dimensionless quantities are kept small. This can 
be established by demanding that the charge of the particle will be much smaller than the 
length scales which characterize the background geometry, i.e., q <^TZ, where TZ = min{7^j}; 
and furthermore we have to assume that e ^ q. To satisfy the last inequality we shall use 
expansion ^ only in the external zone defined to lie outside a worldtube surrounding z{t). 
We denote the interior of this worldtube by S and define its surface S to lie at a spatial 
distance eE^q) from the worldline z{t) such that eE q (e.g., we may take = Nq where 
N is some fixed large number). Later we shall be interested in obtaining the limit q ^ 0. In 
this limit the worldtube radius eE{q) approaches zero, and the perturbation analysis becomes 
valid throughout the entire background spacetime excluding the worldline z{t). 

We now seek a solution to Eq. in the external zone. For this purpose we substitute the 
perturbations series ^ into Eq. ^ and obtain the following equations for the perturbations 

□0(^^ = ,x^S (6) 
□0(2) = v^0WV^0« ,x^S. (7) 

Notice that the form of these equations is independent of q. As they stand Eqs. ()()P7|1 contain 
no information about the physical properties of the charge which induces the scalar field. To 
obtain a unique solution to each of these equations, we must provide additional information 
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about this source. Mathematically, this can be implemented by imposing a set of boundary 
conditions on S which encode the physical properties of the source inside the worldtube at 
each order. 

Before we discuss these boundary conditions in detail, we briefly discuss their physical 
meaning, and provide a rough description of how we intend to impose them mathematically. 
First, let us discuss the first-order equation To formulate the desired boundary condi- 
tions for this equation, we shall study the form of its general solution in a local neighborhood 
of S within the external-zone. Since g -C 7^, there exists a local neighborhood of S in which 
e is much smaller than TZ and at the same time much larger than q; this region is called the 
buffer zone (e.g., in the buffer zone e can be of order \/qTV). By virtue of the smallness of 
in the buffer zone, we may expand the general solution to Eq. (jH)), in powers of 7^~\ 
where the leading order in this expansion is simply the flat spacetime "sum over multipoles" 
solution. By assuming that the source of the scalar field is a structureless point particle, 
we find that its corresponding sum over multipoles can have no higher multipoles than a 
monopole. The coefficient of this monopole term can later be determined from the value of 
the observed charge of the particle. Below we show that by imposing a monopole solution on 
the boundary S we obtain a unique causal solution to Eq. (jH)) in the limit g ^ 0. Next we 
consider the second-order equation ((7j). A general solution of this inhomogeneous equation 
can be written as a sum of a particular inhomogeneous solution and a general homogeneous 
solution. To obtain a unique solution in the limit g ^ 0, the freedom of adding an arbitrary 
homogeneous solution has to be constrained. For this purpose we shall impose particu- 
lar boundary conditions on the homogeneous solution. Here again the absence of internal 
structure in the point-particle source implies that we should impose a monopole boundary 
condition on S. 

Below we shall provide more precise statements about the boundary conditions on S at 
orders q and g^, and we will show how they give rise to a unique causal solution to each of 
the equations (j(j|7|l . in the limit g ^ 0. 

III. FIRST-ORDER PERTURBATION 

Let us focus our attention on Eq. (jH)). Here we shall be interested in obtaining a unique 
causal solution in the limit g — > 0. 
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First, we shall identify the boundary conditions that have to be imposed on S to obtain a 
unique causal solution to Eq. © in the limit of interest. For this purpose, we shall express 
the solution to Eq. ^ using a Kirchhoff representation. In this way the perturbation 0^^^ (x) 
can be expressed as a certain integral over the worldtube surface S. Assuming that the only 
source of the perturbation (p^^^ is the particle itself (i.e. no incoming waves), and assuming 
that 0'-^^ decays sufficiently fast at spatial infinity; we find that we may express 0'-^^ as 

(t)^^\x) = -— [ [G(x|x')V°>(^n^') - 0^'na;')V"'G(x|a;')l . (8) 
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Here G{x\x') denotes a Green's function (which is assumed to fall sufficiently fast into the 
past and the future), the surface integral is evaluated on the worldtube S, and dT,'^ denotes 
the outward directed three-surface element on S. On physical grounds we shall require that 
(f)^^^ will be a causal solution of Eq. To satisfy this requirement it is necessary to impose 
G{x\x') = G^'^^^lxlx'), where (j^^*(x|a;') denotes the retarded scalar Green's function. 

Consider now the general solution 0*^^-' in the buffer zone. For the special case of a fiat 
background spacetime, this general solution reduces to a standard sum over multipoles of 
a static particle. Here the weight of each multipole corresponds to a particular moment 
of the particle's charge distribution. In curved spacetime, it is convenient to describe the 
background metric in a local neighborhood of the worldline which includes the buffer zone, 
using a locally-flat coordinate system based on the worldline. Here we shall employ the Fermi 
normal coordinates. In these coordinates the particle's geodesic trajectory is static. We can 
now solve Eq. (jH)) iteratively. This can be implemented by substituting g^'^ = rj^^ + 5(7^"^ and 
0(^) = 0^^+50^^-* into Eq. (jHl). Here rj^^^ is flat spacetime metric and 0^"* is the corresponding 
general solution in flat spacetime. Notice that in the buffer-zone both 5gf"' and 50(i) are 
small by virtue of the smallness ofe/TZ and q/TZ. We assume that the internal structure of the 
point particle is unaffected by the background curvature. In this case one flnds that S(f){eE) 
is always much smaller than (f)p\eE) as eE 0, and in fact lim^^^o 6(j){eE) /4>e\^e) = 0. 
We can formally express the desired general solution on a family of hypersurfaces r = const, 
generated by spacelike geodesies with tangent vectors which are normal to worldline at z{t), 
as 

oo oo 

0(1) (r, e, = + E • (9) 

m=l m=0 

Here the coefficients {a_m, bm} are independent of e, and Q'^ denotes a unit vector which is 
normal to the worldline and tangent to the geodesic which connects z{t) to x on a particular 



hypersurface. 

By substituting the general solution Q on S into Eq. (jSJ, and recalling that dJ2a scales 
like e^;^, we find that the contribution to the integral which comes from the second sum in 
expansion Q vanishes in the limit q ^ (which implies e£;(g) 0), whereas the first sum 
in Eq. Q which diverges as x -^(t), gives rise to a non- vanishing contribution to (p^^^ 
Note that in the limit q ^ Eq. © is valid throughout the entire spacetime excluding 
z{t), reading 



We conclude that it is sufficient to specify divergent boundary conditions to obtain a unique 
causal solution to Eq. (fTUj) . 

To determine this solution we have to specify the coefficients {a_m} of the divergent 
terms in Eq. Q. First, let us consider the coefficients {a_m} with m > 1. Recall that we 
assumed that the particle has no internal structure. This corresponds to a monopole solution 
in fiat spacetime (i.e., (f)p^ = a_ie^^). Since the spacetime curvature induces corrections of 
a smaller magnitude (as e ^ 0), we find that all the coefficients {a_m} with m > 1 must 
vanish, and the coefficient a_i can have no angular dependence. 

To determine a_i we consider the expansion of the full nonlinear field in powers of 
e in the buffer zone, and identify the coefficient in front of the 1/e monopole term in this 
expansion to be the particle's observed charge. In general this 1/e term may depend on all 
the linear perturbations For simplicity we shall assume that the observed charge is 

time-independent and is identical to q. This implies that a_i = 1 and that the observed 
charge is completely determined by the first-order solution, without any contributions from 
higher-order perturbations. 

Fortunately there is a well-known retarded solution to a linear ho mog eneous wave equa- 
tion (fTn|) with 1/e (monopole) singularity near the particle (see e.g., [id]), reading 



Since our divergent boundary conditions yield a unique causal solution, it must coincide 
with Eq. (ITTD . 

We comment that the same expression pij) can be explicitly derived by substituting the 
expansion ^ with the desired coefficients into Eq. (jHJ, and only then taking the limit 
eg ^ 0. A similar method was used by D'Eath to obtain a solution to the first-order metric 



□0(1) =0 ,x^z{t). 



(10) 




(11) 
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perturbations which are produced by a Schwarzschild black hole P|. In this solution h)^^ in 
the Lorenz gauge takes a form which is very similar to Eq. (fTT|) . 

IV. SECOND-ORDER PERTURBATION 

We now focus our attention on Eq. ((7j). In the limit g — > this equation and its domain 
of validity take the following form 

□0(2) = V^0(^V>(^) , X i z{t) . (12) 

Notice that the form of the source term in this equation is similar to the form of the source 
term in the second-order gravitational perturbations equations there the source term 
contains similar terms which are quadratic in V ah')!}. 

Naively one may try constructing the standard retarded solution to Eq. (|T2|) by integrat- 
ing its source term with the retarded Green's function, reading 

Here the integral extends over the entire spacetime. The difficulty with this naive approach 
is that the integral in Eq. (fT!?|l generically diverges for every point x in spacetime. The origin 
of this divergency is the strong singularity of the source term V^^*-^-* V^^^^-* near the particle. 
Recall that (p^^^ diverges like 1/e, and therefore the second-order source term diverges like 
e~^. Recalling that ^f^^df^x' scales like we find that the expression inside the integral in 
Eq. ()13p scales like e"^^ and therefore this integral diverges. 

Since the retarded solution (fT^ diverges at every point in spacetime it is not a legitimate 
solution to Eq. (fT^ . This difficulty might lead one to believe that Eq. (fT^ has no well 
defined (causal) solutions. Here we show that contrary to this belief there are in fact different 
causal solutions to equation (fT^ which are well defined, and have the desired physical 
properties. Notice that Eq. (fT^ has the following particular solution 

(14) 

this can be easily verified by substituting i/) into Eq. (jl2j) . On physical grounds we have 
to make sure that ip is a causal solution. By a causal solution we mean: a solution to a 
field equation that at every point x is completely determined by the given sources (and/or 
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boundary conditions) of this field equation that are inside and on the past hghtcone asso- 
ciated with point x. Evidently, the retarded solution 0*^^) [given by Eq. (fTT|) ] is a causal 
solution of Eq. (fTUI) . Since ip is completely determined by we find that ip is a. particular 
causal solution to Eq. (fT^ . 

We now discuss the general solution to Eq. (|12p. Consider a different new solution to 

/TTT^ ('2') f2) 

Eq. ()12|) . denoted with 0]y . Defining the difference between 0]^ and our original particular 
solution Tp to be 

we find that (p^^ satisfies the following semi- homogeneous equation (i.e., a homogeneous 
equation away from the worldline) 

Ucf)^" = ,x^z{t). (15) 

Clearly, the general solution to Eq. (fT^ that we denote (pQ^ can be written as 0^?^ = ip+ipQ^ , 
where (p^ denotes the general solution to Eq. 

To construct a unique particular solution (pp to Eq. (fT^ we have to constrain the 
freedom associated with the arbitrariness of the semi-homogeneous part of (p)j . For this 
purpose we shall impose additional set of requirements on (p^^ . These requirements can 
then be transformed to a set of requirements on the semi- homogeneous field (pf^-, yielding 
a unique particular solution to Eq. (fTH|) . which we denote 0p^. By adding this particular 
solution to ip we shall obtain the desired particular solution to Eq. (fT^ . (p^p . 

Note that Eq. p5|) is identical to Eq. (fTUI) . and therefore it can be solved with a similar 
method. Following our previous analysis of the first-order perturbation, we assume that the 
only source of (homogeneous) waves in the fields cp^"^ (and therefore in 0p^) is the particle 
itself. Assuming that decays sufficiently fast at spatial infinity we can express 
using the Kirchhoff representation (jHl). Following the first-order analysis we find that a 
unique causal solution 0p^ is completely determined by its divergent boundary conditions 
as X ^ z{t). Similar to Eq. Q, we formally express (pf>^ in the vicinity of z{t) as 

oo oo 

0^^(r, e, ^]^) = ^ a^le-- + h^^e^ . (16) 

m=l m=0 

Here the divergent boundary conditions for the semi-homogeneous solution (pf,^ can be 
imposed by specifying the coefficients {a^^} of ^^e divergent terms in Eq. (fTBj). 
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We now study the singular properties of in the buffer zone, and from them we shall 
obtain the required coefficients {a^m}- For this purpose we consider g to be a small (finite) 
quantity. In the buffer- zone the first-order scalar perturbation (p^^^ (x) can be expanded on a 
spacelike hypersurface that is generated by geodesies with tangent vectors which are normal 
to worldline at point z{t^), which gives 

,#><'>W=£-'+,#><'>''|j(T.)]+0(e). (17) 

Here. ^<>.« . a regular field wh.ch solves the homogeneous sealar wave equatiou. geuerically 

this field scales like (see jlTl] for the definition and properties of this field). Note that in 
Eq. (fTTj) the field 0^^^^ is evaluated at z(tx) on the worldline. By substituting Eq. (fT7|) into 
Eq. ()14|) we obtain the following expansion 

0g) =^P + <Pf' = \e-' + 0(^)^[^(r.)]e-i + + 0(e°) . (18) 

Substituting this expansion together with Eqs. ()16I17|) into Eq. (0) gives 

= 1^ + [g + g^(0(^)^[^(r,.)] + a^f )] + ^ a^^e"™ + 0(e°) + 0{q') . (19) 

m=2 

As in the first-order analysis, the assumption that the particle has no internal structure 
implies that the coefficients with m > 1 must vanish. We now discuss the coefficient 
of the monopole . Note that the second-term in Eq. (fTU|) has the form of a monopole 
term. We therefore identify the term inside the large square brackets with the observed 
scalar charge of the particle^ (up to second-order). Note that the term 0*-^'*'^[-2(Ta;)] is time 
dependent, and might give rise to a time- dependent observed charge. Since we assumed that 
the particle's observed charge is time-independent and equals g, we find that 

a!f = -0««[.(r.)]. 

Having specified all the coefficients {a^^} of the divergent terms in Eq. (fT^ . the particular 
solution is now uniquely determined. This solution takes the form 



e(x) 



/oo 
0«^[^(r)]G^^*[x|z(r)]c/r, (20) 
-oo 



^ We comment that unlike most linear field theories, the nonlinear field theory which is considered here 
does not have a unique definition for the particle's observed charge. 
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as can be verified by an explicit calculation of the {a^^} coefficients of this solution. We 
conclude that the desired particular solution to Eq. (fT^ is given by 

0g)=^ + 0|^ = -[0«f _ j 0«^[z(r)]G^'=*[x|^(r)]dr, (21) 

where 0^^-* is given by Eq. (|TT| . 

The derivation of the solution 1)211) was based on the following steps: (i) Considering 
the hmit g of the second-order perturbation equation 0. (ii) Guessing a particular 
inhomogeneous second-order solution ip. (iii) Using the physical properties of the particle to 
construct divergent boundary conditions as x — > z{t) for the semi- homogeneous field (p^^ . 
(iv) Solving the semi-homogeneous equation (fT3j) with its divergent boundary conditions 
and obtaining 4>p^. 

Finally we point out that the analysis given here can be generalized to handle the more 
problematic gravitational equation (j21). A detailed account on the regularization of the 
second-order gravitational perturbations will be published elsewhere. 
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